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Abstract 

We present a generator for single top quark production via flavour-changing 
neutral currents. The MEtop event generator allows for Next-to-Leading-Order 
direct top production pp — » t and Leading-Order production of several other single 
top processes. A few packages with definite sets of dimension six operators are 
available. We discuss how to improve the bounds on the effective operators and 
how well new physics can be probed with each set of independent dimension six 
operators. 
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1 Introduction 



The Large Hadron Collider (LHC) at CERN has resumed operation with its center-of- 
mass energy increased to 8 TeV. The LHC top factory will allow us to scrutinise the 
heaviest of all known quarks with unprecedented precision. Flavour physics is on the top 
of the agenda as one of most interesting research topic that can be address at this collider, 
through the study of flavour changing neutral currents (FCNC) in top quark production 
and decay. In fact, a wide variety of models shows a strong dependence in the measurable 
FCNC quantities: for instance, top quark FCNC branching ratios can vary from extremely 
small in the Standard Model (SM) to measurable values at the LHC in a variety of the 
SM extensions pQ. Therefore, the large number of top quarks being produced provide a 
means to search for deviation from the SM, however small they are. It is clear that the 
simplest way to search for new FCNC physics is to look for the rare top decays as for 
example in t — > where q = u, c is the up-quark or a c-quark. Limits on the BR(t — > q^) 
were set indirectly at LEP [2] and HERA0,|I] and directly at the Tevatron[5] (see [6] for 
references and details). Presently, the best bound on the photon FCNC current is the one 
from HERA while the best experimental bounds on BR(t — > qZ) were obtained at the 
Tevatron [7] and at the LHC [H [9] . Finally, the best bound on the strong FCNC current 
tqg was recently obtained in direct top production at the LHC [TO] . 

Our goal is to provide a tool to measure FCNC related to the top quark at the 
production level. More evolved top FCNC searches can only be performed if a dedicated 
generator for top FCNC studies is available. So far only the PROTOS generator [11] 
and the TopReX generator [12] can be used to generate events for direct top production 
pp(gu + gc) — » t studies as well as the top FCNC decays. Therefore, we considered that it 
was both necessary and timely to make available a generator for top FCNC physics that 
included a larger set of FCNC operators together with a larger set of physical processes 
at the production level. This is the purpose of the MEtop event generator. 

MEtop main process is direct top production, pp(gu + gc) — > t, but calculated at Next- 
to- Leading- Order (NLO). NLO direct top was implemented by adopting an effective NLO 
approximation as described in [13J. Besides direct top, MEtop can be used to generate 
events at LO for all FCNC processes with a top and a gluon or any quark other then the 
top in the final state. We plan to include other processes like pp —> tV, with V = 7, Z,W 
in the near future. From the theoretical point of view we will adopt the effective operator 
formalism as described in [H]. We use a set of dimension six effective operators always 
involving at least one top quark. The set of operators is classified in three different groups: 
strong, electroweak and four-fermion (4F) operators. MEtop comes with several choices 
of packages where different sets of effective operators are available. 

The plan of the paper is as follows. The next section describes the complete set of 
FCNC operators needed for top physics. In section 3 we provide a detailed description 
of the physical processes available in MEtop. Section 4 is devoted to the implementation 
in MEtop of the effective NLO approximation for direct top production. In section 5 we 
compare the contributions of the different classes of operators to single top production. 
Our conclusions are drawn in section 6. There are also three appendixes dealing with 
more technical issues. 
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2 The FCNC operators for top physics 



The effective operator formalism assumes that some general theory which has the SM 
as its low energy limit can be written as a series in A with operators obeying the SM 
symmetries, 

C = C™ + + + O (±) , (1) 

where £ 5M is the SM la grangian and £ ( - 5 ^ and contain all the dimension five and six 
operators respectively. This formalism provides a model-independent parametrisation of 
physics beyond the SM. This lagrangian contains only SM fields and therefore any new 
particles and any new interactions are hidden in the effective operators. The term 
vanishes if baryon and lepton number conservation is imposed. 

The complete set of dimension six operators is quite vast. In order to simplify the 
discussion we classify the operators in three categories: strong FCNC operators [151 EE] , 
the ones generating a vertex of the form tug, where g is a gluon and u is an up-quark; 
electroweak FCNC operators [TTJ, dH] , which are the ones giving rise to a vertex with one 
top quark, an up-quark and one electroweak gauge boson, and finally four fermion (4F) 
operators which are Fermi interactions with one top quark and three other quarks. 

When writing all allowed dimension six operators obeying the required symmetries of 
the SM lagrangian, one readily understands that not all operators are independent [H]. 
They are related by the equations of motion and also by Fierz transformations. Therefore, 
the total number of operators can be reduced to a minimum set of independent operators. 
Moreover, this set can be further reduced when only specific processes are studied, like in 
our case, where all operators have at least one top quark in the interaction. A minimal 
set of operators for top quark physics was discussed in [T5J [HI (El 120] and here we will 
just present this minimal set according to our classification. 

We will start by considering the non 4F operators. As previously stated, we divide 
these operators in two classes: strong FCNC operators, when the gluonic tensor is in- 
volved, and electroweak FCNC operators when electroweak gauge bosons are present in 
the interaction. We assume that O y and O^ 1 are independent operators and the hermitian 
conjugate of all the operators are included in the final lagrangian. 

Following the notation of [2], the independent operators contributing to the strong 
FCNC vertices can be written as 

0% = <? L \ a o» v u^G^ , (2) 

where G° u is the gluonic field tensor, u l R stands for a right-handed quark singlet and q l L 
represents the left-handed quark doublet. FCNC occurs because one of the indices is 
always equal to 3 while the other is either 1 or 2, that is, there is always one (and one 
only) top-quark present in the operator; the remaining fermion field in the interaction 
is either a u or a c-quark. These operators will give rise to the FCNC vertices of the 
form gtUi (with m = u, c) and the corresponding hermitian conjugate interaction with 
an independent coefficient. Operator in equation (J2J) also appears in the literature as a 
dimension 5 operator. In that case, the corresponding FCNC lagrangian is written as 

C s = iK u il u \V(f u + h u 75 )tG" + (u c) + h.c. (3) 
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where k u is real, g s is the strong coupling and f u and h u are complex numbers with 
\fu\ 2 + \ h u \ 2 = 1 (see appendix A for a detailed discussion relating the forms of the strong 
FCNC operators). 

In the electroweak sector we now have to look for FCNC vertices of the type VtUi 
(with Ui = u, c and V — Z, 7). The minimal set of operators that give rise to the above 
interactions can be written as 

0% = t ^ v> R B, u , = t rj o>» v? R } Wl , (4) 

0% = i^D^^ui), (5) 
Oft* = * ItfDJ) (<f L Y ql) , = i (0t Tl D,cj>) {q\ >f n <f L ) , (6) 

0% =(tf<t>)(qiuU), (7) 

where B^ u and W* are the U{l)y and SU{2)l field tensors, respectively. As was shown 
in [19], for all the operators in fl5]) and fl6]), and are not independent. In fact, by 
writing the combinations O t+ ^ and O l ~i and using the equations of motions, it can be 
shown that only one of these combinations is independent. This means that the number 
of independent operators in 03) and (EJ) is reduced to three (for each light flavour). The 
above discussion leads us to the conclusion that the minimal number of non-4F operators 
needed to study top FCNC physics is 9 for each light flavour {u and c). 

The equations of motion used to reduce the number of operators relate the operators 
from the strong and electroweak sectors with the 4F operators. The number of indepen- 
dent 4F operators depends on the process considered. Following [2T] we found the minimal 
number of independent 4F operators needed for t plus quark production, pp — > tq, which 
are shown in appendix B. The final lagrangian for the study of single top production via 
FCNC currents can then be written as 

r> ' \ v ( ij /nij I ij /nij , ij /nij . ij /nij . {3,ij) /n(3,ij) 

L, 1<im,99- J ftq — ^2 y^uGip U uG(/> ' a uW<t> U uW<t> ' a uB<f> U uB4> ' a <f>u U 0u ' a 4>q U 4>q 

»J=1)3 

or 
i,j=2,3 



0& w + a u , Ol) + - C Afu + - £ 4/c (8) 

where £ 4 / u and £ 4 / c are described in appendix B. One should note that in order to keep 
a manageable number of 4F operators we only consider initial states with up-quarks in 
the hadron colliders case. There are reasons for considering a reduced set of 4F operators, 
namely the ones that have in the initial state either uu and uu. First, and assuming 
that all 4F coupling constants are of the same order, these initial states provide the 
largest contribution for the cross section. Second, our main goal is to provide a means 
to distinguish between operators by analysing different distributions and this can only be 
done if the number of operators is not too large. The addition of further 4F operators 
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will be done in the future if found necessary. The minimal number of 4F operators in the 
case of FCNC pp — > ti was recently considered in [22J . 

Before ending this section we will briefly discuss the bounds on the coupling constants 
«j. In the effective operator approach, the lagrangian has the SM symmetries. Therefore, 
physics of the top quark is related with B physics. In reference [23] this relation was 
explored in order to constrain the electroweak FCNC operators Q. The most constrained 
operators are obviously the ones built with quark doublets only while the less constrained 
are the ones built with quark singlets only. Consequently, for the first generation, bounds 
on operators with doublets only, och /A 2 and a^^/A 2 , are of the order 0.01 TeV -2 . 
For operators a^^/A 2 and ct^ Br p/ A the bounds are of the order 0.3 TeV -2 while for 
a^^/A 2 and a^^/A 2 we have 1 TeV -2 . Finally, regarding operators with singlets only, 
like a u< f)/A 2 , reference [23J obtained a bound of the order 3 TeV -2 . The bounds for 
operators relating the second and third generation are of the same order of magnitude. 

As stated in the introduction, there are new direct bounds from the LHC that lowered 
the limit on BR{t — > qZ) to 0.34% [9J. A new indirect bound from HERA [3] is also 
available BR(t — > q'j) < 0.5%. Also, a combined study on B physics and Tevatron data 
on top quark production cross section places an indirect bound on the sum of the FCNC 
branching ratios forcing them to be below the percent level [25J. All these new bounds 
do not imply any dramatic changes on the bounds in the electroweak sector. 

Regarding the strong FCNC operators the most recent search is the one from the 
ATLAS collaboration [10] in direct top production at the LHC. The upper limits obtained 
at 95 % CL for the strong couplings are k u /A < 6.9 x 10 -3 TeV -1 and k c /A < 1.6 x 10 -2 
TeV -1 which in turn can be translated into strong branching ratio bounds BR(t — > ug) < 
5.7 x 10 -5 and BR(t eg) < 2.7 x 10 -4 . 

Contrary to the the dimension six FCNC operators from the strong and electroweak 
sector, there are no useful bounds on the four fermion operators involving two top quarks 
and this is even more so if the top is right-handed. Therefore, the LHC can place con- 
straints on these operators. 

3 Physical processes 

MEtop generates events according to the Von Neumann algorithm (see [26j for details). 
The amplitudes for each process were generated with CalcHEP [27J, and the Feynman 
rules for the effective operators were derived with LanHEP [28J. Integrations are per- 
formed using the CUBA library |29| . configured to use VEGAS algorithm [30] Generation 
of events for hadron colliders need the linking with the LHAPDF package [31J. The events 
are written in the standard format of Les Houches event file [32] (.LHE). Whenever pos- 
sible our results were checked via a completely different path. First, the Feynman rules 
were generated by the implementation of the effective operators in UFO [33]. Then, cross 
sections calculation and event generation was performed using MadGraph 5 [34J. We 
always found an excellent agreement with MEtop. 

The following hard processes are already included in MEtop: direct top production 

1 Other analysis based on B physics observables and electroweak precision constraints were also per- 
formed in [24] leading to similar conclusions. 
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at LO and NLO both for pp and for pp colliders, which at the parton level amounts to 
the processes gq — >■ t and gu — > gt where q — u, c; pp{p) — » gt, with all possible parton 
contributions in the initial state taken into account and q is now any quark other then 
the top-quark. The corresponding conjugate processes, with an anti-top in the final state, 
are also included, and the processes can be generated independently. There are several 
available packages in MEtop, each containing a different set of operators. The list of 
operators in each package is presented in appendix C. 

In figure [1] we present direct top production together with the top+gluon processes. 
The two have in common the fact that only the strong FCNC couplings contribute to 
the process. Only one diagram contributes to direct top channel while several diagrams 
with gluon exchange contribute to top+gluon. MEtop allows for generation of events at 
LO and at NLO. As we have already mentioned, currently, FCNC direct top production 
events can be generated at LO with the PROTOS generator. 





9 9 t 

Figure 1: FCNC leading order direct top production and top + gluon production at the 
parton level. Only FCNC strong operators contribute to the process. 



In figure [2] the diagrams for top + quark production are shown. Both qi, q2 and q 
run through all quarks other then the top-quark, that is u, d, c, s, b and the respective 
anti-particles. Again, conjugate processes are also included. In the diagram on the right, 
only strong FCNC operators are present. On the left diagram all operators can take part, 
including the 4F ones, contributing to LO single top production at the parton level. A 

qi q 9 q 




92 t 9 t 

Figure 2: FCNC leading order top + quark production at the parton level. FCNC strong 
and electroweak operators contribute to the process together with 4F operators. 



detailed description of each process will be presented in appendix C. 
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4 NLO approximation to direct top production 

In an event generator, initial (ISR) and final (FSR) state radiation is included trough a 
shower mechanism. The shower mechanism assumes collinear factorization, that is, that 
the real radiation process can be approximated by a branching mechanism, where the first 
QCD radiation is emitted by one of the legs of the Born configuration. In equation [9] we 
present the relation between the transition amplitude for the case of q — > qg splitting 



where Mj is the amplitude and d<&i is the phase space for the zth body processes and 
Pq,q g (z) is the Altarelli-Parisi splitting function. This approximation breaks down in 
the hard Pt region where the matrix element M n+ \ should be used. A factorization 
prescription or matching scheme is then used to merge these two regions in a smooth and 
optimised way (see for example the merging approaches CKKW [35J and MLM [36J). 

In the previous section we have presented in figure Q] the parton level contributions 
to direct top production, together with the parton level contribution to the hard process 
gq —7- gt. The later process contributes to the inclusive direct top production. The 
gq —7- gt process has soft and collinear divergences, and this problem can only be solved 
by including the NLO corrections. Furthermore, the FCNC direct top production cross 
section was calculated in [37] and a considerable enhancement of about 40 % was found 
relative to the leading order cross section for the LHC with yfs = 14 TeV. Therefore it 
is desirable to have a NLO generator for direct top production at the LHC. It was also 
shown in [37] that the NLO QCD corrections vastly reduces the dependence of the total 
cross section in the renormalization and factorization scales which in turn increases the 
confidence in the predictions. 
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Figure 3: Inclusive FCNC direct top production at NLO in QCD. 
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In order to calculate the inclusive NLO FCNC direct top cross section one has to 
consider the Born diagram, the virtual contributions and the real emission diagrams. 
A sketch of the diagrams from each of the above mentioned contribution are shown in 
figure [H It is well known that at NLO, the integration in the full phase-space gives rise 
to infrared divergence from the virtual-born interference part and these divergences are 
cured by the addition of the real emission process. Although easily fixed for a total cross 
section analytical calculation, this problem is not straightforward to deal with at an event 
generator level. In fact, there is no way to extract the infrared infinities as it is usually 
done with dimensional regularization in a theoretical calculation. There are methods 
to deal with these infinities such as Phase Space Slicing [38j EHl HQ] and Subtraction 
Methods [HI 02]. In this work we will adopt an effective NLO approximation [13] to 
simulate direct top events at the NLO level. In this approach a merging scheme between 
2 — > 1 and 2 — >■ 2 events is performed, where each process will separately populate two 
distinct but joint regions of the phase-space. A resolution parameter must then be defined, 
which in the present case is the cut in transverse momentum of the top quark applied to 
the real radiation process. This P T ut will then play a role of a matching variable, p™ atch . 
The phase-space region for small Pt will be described by the 2 — > 1 process and the 
subsequent parton shower (PS) mechanism, whereas the hard Pt region will be described 
by the 2 — > 2 process. One must then just make sure that the transition is done in a 
smooth way. The virtual corrections are included via a K-factor applied to the cross 
section of the 2 — >■ 1 process. We assume this to be a good approximation because the 
kinematics of the Born and Virtual configurations of the direct top process should be very 
similar. The events will then be generated according to the following relation 

a NLO = Ka LO (Pf S < Pr tch ) + VReai(P T > Pr tch ) (10) 

where a to is the tree-level direct top contribution, aR ea i is the real radiation part, K is the 
K-factor and P^ s and p™ atch are the transverse momentum of the first PS emission and 
the integration cut of the real radiation process, respectively. Once the direct top events 
are produced, they will be radiated through a radiator like the one in PYTHIA [33]. In 
order to avoid double counting, the matching must ensure that the first PS emission from 
the 2—7-1 process will not fall within the 2—7-2 configuration phase-space. There are 
two ways of accomplishing it: either by vetoing all radiated 2 — > 1 events that would be 
within the 2 — >■ 2 configuration phase-space or simply by limiting the phase-space region 
of the radiated 2—^1 events to the boundaries defined by the resolution variable. We 
choose to adopt the later. 

In order to follow this approach, one must ensure that the PS mechanism added to the 
generated events from the Born configuration will populate the region with Pt < P T natch i 
which can be assured using a PT-ordered shower [33] , available in both current PYTHIA 
versions 6.4 and 8.1. We therefore assume that the generated events will be showered by 
a PT-ordered mechanism. Therefore we start by calculating the three cross sections from 
equation (fIU|) . with P T ut = p™ atch for the 2—7-2 process. For the (Jnlo cross section we 
have used the expressions from [37J, where the top quark is on-shell. The tree- level direct 
top and top+gluon amplitudes were generated with CalcHEP where the top quark and the 
W decays were included in order to preserve spin correlations. The cross sections are then 
calculated with the Cuba library. Hence, the K factor is calculated "on the fly" for each 
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sub-process. After extracting the K-factor, the events are generated weighed according 
to equation [TOJ The PS starting scale can then be configured to start the branching in 
pmatch f or 2 _ y \ eV ents configuration, which in MEtop is done by preparing the .lhe 
files to be used by PYTHIA. A short remark is in order - in the 2 — > 1 configuration, 
no meaningful FSR from the top quark can be present due to its large mass. Hence, we 
consider a good approximation to take only ISR into account. 

f Z\top after first ISR I No Match 

1 dN 

(GeV" 1 ) 

N dP T 

0.14 t r , 

0.12 || 
0.10 ;!] 
0.08 | '; 
0.06 | \ 
0.04 ! V 




20 40 60 80 100 



p p -> t 

p p -> t+g 



Figure 4: Pt distribution of the top quark for y/s = 7 TeV. The black solid line is for 
direct top production after the first branching in ISR, with starting scale of m t . The blue 
dashed line is for the hard process top+gluon production. 



In figure HJ the black solid line represents the Pt distribution of the top quark in 
direct top production, after the first branching in ISR, with starting scale of rrit- In the 
same figure, the blue dashed line represents the hard process: top+gluon production. As 
described previously, Pt is the kinematical variable chosen to match the two processes 
avoiding double counting in the low Pt region. 
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Figure 5: Pt distribution of top quark after the first ISR branching with a p™ atch of 10 
GeV (left) and 15 GeV (right). 

In figure [5] we present the Pt distribution of the top quark after the first ISR branching 
with a p™ atch of 10 GeV (left) and 15 GeV (right). The natural criterion to determine the 
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value of the Pt matching parameter in the effective NLO approximation is the smoothness 
of the transverse momentum distribution. There are no significant differences when the 
value of p™ atch is varied in the 5 GeV to 20 GeV range. As can be seen in figure [5], there 
is never a completely smooth transition between the two sets of events. This effect should 
be included as part of the systematic uncertainties. This feature was checked for a large 
range of Pt match. After including the full shower (ISR+FSR) and Multiple Interaction 
(MI) we have opted for a value of P T natch of 10 GeV. 



Pt, top isr+fsr+mi I Match = 10 GeV 




Figure 6: Pt (left) and r\ (right) distributions of the top quark at the partonic level after 
the full shower (ISR+FSR) and Multiple Interaction. 

In figure [6] we show the Pt (left) and r\ (right) distributions of the top quark at the 
partonic level after the full shower and MI for p™ atc/l = 10 GeV. The blue dashed line 
represents the real radiation part while the grey dashed line is the direct top fully showered 
but with the pt starting scale at 10 GeV. The solid black line is the final NLO distribution 
which amounts to the sum of the previous two. 

^V.top '/top 
1 dN 1 dN 
(GeV -1 ) 




Figure 7: Comparison of the LO and NLO Pt (left) and r\ (right) distributions of the top 
quark at the partonic level after the full shower (ISR+FSR) and Multiple Interaction. 

In figure [7J we present the LO and NLO Pt (left) and rj (right) distributions of the 
top quark at the partonic level after the full shower and MI. It is clear from figure [7J 
that the Pt and rj distributions of LO direct top production are quite different from the 
corresponding NLO direct top ones. In fact, the distributions show that the use of a 
constant K factor does not correctly describe the behaviour of direct top at NLO. Hence, 
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a new analysis is needed to improve the accuracy of the bounds on the strong coupling 
constants k u and k c . The direct top NLO pt distribution is shifted to low values of px as 
compared to the LO distribution while the rj distributions are shifted to higher values of 
i] as compared to the LO one. 

The actual experimental analysis is performed by looking at the distributions of the 
final state particles. Therefore, in figure M we present the comparison between LO and 
NLO Pt (left) and rj (right) distributions of the lepton from t — > bW — > blv at the 
partonic level after the full shower and MI. Again, it is clear that the level of improvement 
by considering the NLO distributions heavily depends on the particular analysis being 
performed. 



T , lepton 



'/lepton 



1 dN 

(GeV -1 ) 

NLO dPj- 




Figure 8: Comparison of the LO and NLO Pt (left) and rj (right) distributions of the 
lepton from t — > bW — > blv at the partonic level after the full shower (ISR+FSR) and 
Multiple Interaction. 



Finally, in figure [9] we compare the LO and NLO Pt (left) and 77 (right) distributions 
of the b-quark coming from t — > bW — > blv at the partonic level after the full shower and 
MI. 




Figure 9: Comparison of the LO and NLO Pt (left) and rj (right) distributions of the 
b-quark from t — > bW — > blv at the partonic level after the full shower (ISR+FSR) and 
Multiple Interaction. 



We have just described how we generate a sample of inclusive direct top production at 
NLO. However, if the goal is to set a limit on the strong FCNC coupling, one needs to add 
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the events generated in the process pp — > t + jet composed by the parton level processes 
gg —> t + q and qiq 2 — > t + q, where q, qi, q 2 are quarks other then the top and hermitian 
conjugate reactions are included. Although the main contribution to top + jet production 
comes from the reaction where the jet is a gluon, all processes where the strong FCNC 
operator intervenes should be taken into account in all analysis. As these processes also 
suffer from infrared and collinear divergences, we have decided to avoid them by using a 
similar cut to that of p™ atch ( that is, pt > 10 GeV. The complete NLO QCD corrections 
to the FCNC process of top+jet production were presented in [15]. The corrections can 
increase the cross section by 10 % to 30 % at the LHC@14TeV. 

When generating the top + quark subprocesses we have to decide what is considered 
as signal in our analysis, pp — > tq has three different classes of subprocesses: the ones 
which are exclusive to the Standard Model, like ud — > tb, the ones that are originated 
exclusively via FCNC interactions, e.g. uu — > tu, and the ones where interference occurs, 
like ub — > tb. We define as FCNC signal the contributions from the two later classes 
of subprocesses. For the pure FCNC processes this poses no problems. However, for 
the interference terms this procedure leads to the inclusion of a small portion of events 
that will also be counted as background. However, choosing the effective strong coupling 
constants as n u = 0.01 (A = 1 TeV) , f u = l/y/2 and h u = l/y/2 and for a CM energy of 7 
TeV, the pure FCNC cross section is 8.718 pb, the interference term is 1.205 pb while the 
SM contribution amounts to only 0.018 pb. Hence, the SM contributions can be safely 
neglected. 



f*T,top Vtop 




NLO 
NLO+'t+q' 



Figure 10: Pt (left) and rj (right) distributions of the top quark for NLO direct top (solid 
line) and NLO direct top plus pp -» tq with p™ atch = 10 GeV and jet p T > 10 GeV. 

In figure [10] we show the Pt and rj distributions for the direct top at NLO summed 
with pp — > tq for a p™ atch = 10 GeV and the jet pt > 10 GeV. It is clear that the shape 
of the distributions do not change much with the inclusion of the pp — > tq process but still 
the pp tq process gives a contribution of the order of 10 % to the total cross section of 
the inclusive top production at the LHC at = 7 TeV. 



5 Single top beyond the strong FCNC operators 

In the previous sections we have discussed NLO direct top and t + q production when 
only the strong FCNC operator is considered. We note that the leading order contribution 
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to direct top does not receive contributions from other operators. Therefore, the NLO 
calculation is again performed with only the strong FCNC operator. As long as no excess 
is found at the LHC in the single top channel, the procedure described in the previous 
section gives us the best possible bound on the anomalous strong FCNC coupling when 
all other operators are discarded. 

In the hard Pt region, the process pp — > t + jet gets contributions from the complete 
set of independent operators. As these operators are independent from each other (and 
therefore so are the respective coupling constants) the interference terms between strong 
and electroweak or 4F could be sizeable. If an excess is found in the single top channel, 
one has to take into account all possible contributions from the remaining operators. A 
thorough analysis of the distributions of each individual operator will help us understand 
which operators could be important for a given experimental analysis. Moreover, even if 
an excess is not seen in the single top channel, dedicated analysis could most probably 
help constraining definite sets of operators. 




Figure 11: Pt (left) and 77 (right) distributions of the top quark when only the strong 
operator is turned on with P T ut = p™ atc/l = 10 GeV. Process considered is direct top at 
NLO plus pp — > tq for a/s = 7 TeV and three values of k u with A = 1 TeV. 



We start by considering the strong operator. When all other operators are turned off 
the Pt and rj distributions have a very mild dependence on the strong coupling constant k 
(k u to be more precise, and we have set k c — 0). This is shown in figure ITT1 where the Pt 
(left) and 77 (right) distributions of the top quark are shown for three values of k, 0.001, 
0.01 and 0.1 and A = 1 TeV. The process is direct top NLO plus pp — > tq for y/s = 7 TeV 
and P T ut = P T natch = 10 GeV. It is clear that the shape of the distributions does not vary 
much making it possible to perform the analysis for one constant and then to extract a 
bound on the strong operator. 

We now move to the study of the electroweak operators. We first consider only one 
operator O u w<p turned on. In figure IT2l we present the Pt (left) and i] (right) distributions 
of the top quark for three values of a u w4>, 0.01 , 0.1 and 1 and A = 1 TeV. As a u w4> — > 
we recover the pure SM contribution of electroweak origin. The SM cross section for this 
process and for 7 TeV is a = 0.019 pb while the total cross section for a u wcj, = 0.01, 0.1 
and 1 are a = 0.0020 pb, 0.148 pb and 12.4 pb respectively. Therefore the different shapes 
of the Pt and rj distributions are due to the interference with the SM contribution. When 
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Figure 12: Pt (left) and rj (right) distributions of the top quark when just one electroweak 
operator, O u vy</>, is turned on. The process is pp — > tq for y/s = 7 TeV and P T ut = 10 
GeV. 

ocuW(j> — 1; the total cross section is almost 100 times larger than its pure SM counterpart. 
Therefore, this value shows how the distribution behaves when the SM contribution is 
negligible. 

This kind of behaviour can occur for any operator on the list, provided that the 
coupling constants are such that SM and FCNC cross sections are of the same order 
of magnitude. Any deviation relative to the SM showing up in the distributions could 
mean an interference with one or more operators. Understanding the different distribution 
requires dedicated studies with no assurance however that the responsible operators could 
be identified. One should emphasise that a thorough study of the Pt and rj distributions 
of the top quark could help identifying classes of FCNC operators. A similar discussion 
applies to the 4F operators case. 

One should also note that pp — > tq does not include direct top, because the strong 
operator is turned off. Contrary to strong operator scenario, in this case the distributions 
change with the value of the electroweak constant. Therefore, any bound based on the 
production process has to take into account that different coupling constants can lead to 
different distributions. 

Having studied the distributions of a definite operator representative of each class, we 
will now perform a comparison between classes. In figure [13] we present the Pt (left) and 
rj (right) distributions of the top quark when just one operator is taken non-zero at a time. 
We compare the distributions of the strong FCNC operator with one electroweak operator 
(with coupling constant at u w<t>), an d one 4F operator, (u7 M 7^ u) («7 M 7l t). It is clear from 
the figure that the distributions can be quite different and therefore distinguishable to 
some extent. The ability to distinguish the different operators depends heavily on the 
relative values of the coupling constants. If an excess in single top production is seen we 
can try to understand its origin by looking at all possible distributions. However, this will 
always be a hard task because different operators give similar distributions and therefore 
only very particular scenarios can be probed. 
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Figure 13: Pt (left) and 77 (right) distributions of the top quark when just one operator 
is taken non-zero at a time. We compare the distributions of the strong FCNC operator 
with one electroweak, O u w<f>, and one 4F operator. The process is pp — > tq for y/s = 7 
TeV and P T ut = 10 GeV. 

6 Conclusions 

We have presented a new generator for the study of FCNC top interactions. The generator 
MEtop comes with different packages, each with a subset of a complete set of dimension 
six operators. At the moment MEtop can generate events for direct top and for top plus 
jet production, where the jet can be any quark other then the top or it can be a gluon. 

The direct top production process is implemented at NLO using an effective NLO 
approximation. Also, the inclusive contribution to direct top coming from pp — > tq can 
be included in the event generation. We have shown that the top Pt and rj distributions 
show clear differences when the events are generated at LO or at NLO. Therefore, the use 
of a constant K-factor does not provide an accurate description of direct top production at 
NLO. We conclude that a new experimental analysis is in order to improve the constraints 
on the strong FCNC coupling constants. The inclusion of the inclusive process pp — > tq will 
further improve this bound. We note that a detailed study of the Pt and 77 distributions 
of the top quark could help identifying classes of FCNC operators . 

At LO, the contributions stemming from the different operators can be compared in 
the single top production process. In particular, 4F operators can be for the first time 
constrained at hadron colliders. Constraining the 4F operators can help us understand 
their role in the asymmetry measured at the Tevatron. 

The bounds on BR(t — > u(c)Z) and BR{t — > u(c)j) are obtained in the process 
pp — > ti where one of the top-quarks decays as t — > bW while the other decays as 
t — > u(c)Z or t — > u(c) , ~f. This means that all electroweak FCNC couplings always 
appear in the same combination. With MEtop we are able to look for distribution that 
isolates each electroweak FCNC operator. This way more detailed information can be 
obtained about each operator. 

New final states with FCNC contributions, like for instance pp — > tW [UJ , are to be 
included in the next version of MEtop. 
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A Effective couplings translation 

The representation of the effective coupling constants for each operator is arbitrary. In 
this appendix we will relate the two most common representations of the strong FCNC 
coupling constants appearing in the literature. The translation is simple and our goal 
is to clarify the relation between experimental bounds (and theoretical bounds as well) 
present in the literature. When the effective strong operator comes from a dimension six 
effective lagrangian it is usually represented in the form 

4 6 = ^Ql^^nUG^ +h.c, (11) 

where A is the scale of new physics. This operator can be written as a dimension five-like 
operator when the scalar field, is replaced by (v/y/2 0), resulting in 

= ^°%u\ \ a u j R G a » u + h.c. . (12) 

There are in principle four complex constant involved, an and au with i = u, c in a total 
of four degrees of freedom. On the other hand, several authors adopt to write the same 
strong operator as a dimension five effective operator. In this case it is usually written as 

C V = -ds^u 1 X a (/ fj + ih ijl5 ) v? G a ^ + h.c, (13) 

where gs is the strong coupling constant, «y is taken as real and positive, /y and 
are complex and \fij\ 2 + \hij\ 2 = 1. In both cases Tr[\ a \ b ] = S a b/2, and the vacuum 
expectation value is v = 246 GeV. Note that in reference [H] (as in other references) the 
vacuum expectation value is defined as v — 246/ y/2 GeV. The constant real and 

only two constants are needed, one for each light flavour; the same is true for the complex 
constants f\j and hy. Therefore we can use just the light quark index to represent those 
coupling constants 

Cf = -gsju* X a a» v {fi + ih ilb ) v? G a ^ + h.c, (14) 
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and in this case the index j refers to the top quark. It is now straightforward to find the 
relation between the two set of operators 

an = -V2g s — K i (f i + ih i ) 
v 

a ti = -V2g s - Ki (f* + ih*) (15) 
v 

with % — u,c. In most cases all constants are taken as real. This means that equation [14] 
can be written as 

Cf = -g s ju* X a o»» (fi + h a ,) v? G^ v + h.c, (16) 

and consequently 

a it = -V2g s — Ki (fi + hi) 

v 

a u = -V2g s — K i (fi-h i ) (17) 
with % = u,c, all constants are now real and \fi\ 2 + \hi\ 2 = 1. 

B The complete dimension six lagrangian for single 
top production 

As previously discussed, the number of effective dimension six operators is huge. There- 
fore, no meaningful analysis is possible when all operators are considered simultaneously. 
Any subset of operators, however small, does not simplify much the task of obtaining 
information about each individual operator. In order to perform any relevant study in- 
volving all different types of 4F operators, we have built a 4F subset to be used as a 
basis for our study. Since we are dealing with hadron colliders, our first simplification 
is to consider only processes initiated by u-quarks. This is equivalent to say that the 
coupling constants in the 4F sector are all of the same order. Using references [T41 1201 [21] 
we extract nine 4F effective operators that could contribute to single top production at 
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hadron colliders. This set can be written as 



/nijkl _ 

qq 


1 

^{QLii^qLj) (qLkiuqu), 


(18) 


/nijkl _ 

qq' 


^ {qLia Y qijb) {qikb la qua), 


(19) 


/nijkl _ 
^ uu 




(20) 


/nijkl 
°ud = 
/nijkl 
U ud> — 

(nijkl 

^ qu 

/nijkl _ 

qu' 
/nijkl 
W qd — 
/nijkl 
U qd> — 


(u Ri Y u Rj ) (d Rk 7 M d m ), 

(uRiaY u Rjb) (d~RkbJv,dma), 

{quURj) {uRkqu), 

((iLiaURjb) (uRkbqLla), 

{qu dRj) (d Rk q L i), 
{quadRjb) (dRkbqiia), 


(21) 
(22) 
(23) 
(24) 
(25) 

(26) 
(27) 



where i, j, k, I — 1, 2, 3 are flavour indices and the sub-indices a and b indicate the contrac- 
tion of color indices whenever this pairing is different from the one in spinor contraction. 
In [21] this set of operators was simplified to the one presented in table [TJ There is a total 



|«f + ^ 6 ){ULkYu L j){uul,t L ) 


-\ a m' 3 (uLkYu Lj )(u R a^t R ) 


~ l a w' d (uRkYURj ) {Uul^L ) 


\a k J^(uR k ^URj)(u R a u tR) 


-^Oif* 3 (ULkaYULjb) {uRib1^tR a ) 


-^qu\uR ka ^UR j b){u Lia ^^t La ) 



Table 1: 4F operators contributing to single top production with parton level processes 
of the type u(u)u{u) — > t(t)u(ccu). 

of 24 different combinations coming from setting two of the indices i, j, k equal to 1 while 
the remaining one is set to 1 or 2. By forcing the initial state to be composed of u-quarks 
only, we can further reduce the number of operators to 12 (we only allow for one "FCNC- 
current", (uTt) or (cTt), where T stands for a generic Lorentz structure). Therefore the 
final 12 independent 4F operators are the ones obtained by setting k = j = 1 and i = 1,2. 
Because the operators had to be rearranged to allow for the implementation in LanHEP 
(see discussion below), the final lagrangian is then written as 

f ^ \ * ( ij /nij i ij /nij i ij /nij i ij /nij i (3,ij) /n(3,ij) 

L, qqm,99^q. ~ ^2 \ a uG^uG<l> ' a uW<j> U uWcj> + a uB<f> U uB<j> + a <f>u U cj>u + a <f>q U (j>q 

i,i=l,3 
or 

+4f Off + a u , 0%) + i C Afu + i £ 4/c (28) 
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where £ 4 / n is the 4F lagrangian for anomalous top-up coupling 
Ufu = ^« 13 + a 1 ^){u L ^u L ){u Ll ^t L ) 

-\Ku> 3 + l^ q u 3 )(u R Yu R )(u Llll t L ) 

+ \a^u\u R Yu R ){u R -iJ R ) 
-\a\f\u L ^\ a u L ){u Rl ,\ a t R ) 

-\a^{u R ^\ a u R ){u Ll ,\ a t L ) (29) 
and C^fc is the 4F lagrangian for anomalous top-charm coupling 

-^Hl* 1 + \a™ l )(u Ll »u L )(c Rl J R ) 

-\a^\u L r\ a u L ){c Rlll X a t R ) 
-\af^(u R Y\ a u R )(c Ll ,\ a t L ) 
+ ^« 13 + aH? 3 )(c L YX a t L )(u L YX a u L ) 

+ \all 13 (c R Y\ a t R )(u R Y\ a u R ) 

I i 

+ (^ 2 q l} 1 + a 2 q l n )(c L t R )(u R u L ) + -a 2 q l}\c L \ a t R )(u R \ a u L ) 

+ (^? 3 + a^){c R t L ){u L u R )) + \a)$\c R \ a t L ){u L \ a u R ). (30) 

The operators with the Gell-Mann matrices originate from re-writing the ones where the 
quark colours indices were explicitly summed. The inclusion of 4F operators in MEtop 
was done by implementing the 4F effective lagrangian in LanHEP. All 4F operators in 
table [T] have four coloured particles converging in one point which is a type of interaction 
LanHEP is not able to handle automatically due to the complex color flow. Therefore, 
we had to implement these operators using an auxiliary field mechanism [46], where the 
4-color vertex is replaced by 3-color vertices that when combined in s,t and u channels, 
will reconstruct the 4-fermion interaction. These 3-color vertices are implemented by 
introducing the interaction terms in the initial lagrangian together with a unit mass field 
with a point-like propagator. An example of how a lagrangian is written is shown in 
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equation EH with a vectorial auxiliary field 

uf = {^i^mhM (^7^D^ + + 5V ( 3l ) 

where V'l is a left-handed spinor and is a spin 1 field that does not propagate. 

C Using MEtop 
C.l Installation 

MEtop is written in C and python and it generates events following the Les Houches 
Accord format. It can therefore be easily interfaced with PYTHIA or Herwig. In order to 
compile it, you need a C compiler and python version 2.6 or later. To run the package 
you must additionally install 

• Cuba Library version 3.0 

• LHAPDF version 5.8.6 

• Numpy version 1.3.0 

The Cuba and LHAPDF library must be available through the library environment vari- 
able (for example). 

To install MEtop you just have to execute "make" in the main directory. 

C.2 The generator 
C.2.1 param.dat 

In MEtop all parameters are set in one file: "param.dat". Table [2] summarizes the defini- 
tion of each parameter. 

C.2. 2 Physical processes 

In addition to the parameters defined in table |2] there are two more flags in "param.dat" 
file: "cs" and "Process". The first one dictates whether or not to calculate the cross 
sections and/or to generate events. The second sets which physical process should be 
taken into account. If "cs" is set to 0, the cross sections for all sub-processes defined by 
the "Process" flag will be calculated and no generation will be performed. The result 
will be stored in the CS folder, in a csX.txt file, where X can be "Dtop","Gtop" and 
"Lqtop". If "cs" is set to 1, only the event generation will be performed. In this case 
events are produced according to the calculated cross sections. After generation, the LHE 
files will be stored in the Events folder together with a file "runinfo.txt" which stores all 
information related to the event generation. 

2 There is one file written in Fortran and therefore you also need a Fortran compiler. 
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Mx 


Particle's masses (x=u,d,c,s,b,top,e,/i,r,W,Z,H) 


wx 


r articte s Wicltns ( k x= W,top,Zj,riJ 


sx 


Values for CKM matrix elements (x=12,23,13) 


b W 


singly ( t>w is the Weinberg angle) 


EE 


Electromagnetic coupling constant 


cox 


couplings of the x operator (x=l,2,...,9) 


C 1 

ix, nx 


Chirahty parameters from operators col and co2 


Q 


T7 1 J- * J- * 1 

1 actonzation scale 


miuR 


Renormalization scale for Direct top at NLO 


T 

L 


Energy scale 




Centre of mass Energy 


PTmatch 


PT for matching 


PTmin 


Cut in PT for LO 2 — > 2 processes 


INrL/vnts 


Number of events to generate 


pat 


PDF name according to LHAPDF 




Tvoe of collider - 1 for t>t> and -1 for t>t> 


DecMod 


Turn on/off W decay modes 


SpCorr 


Turn on/off Spin Correlations 


ttbar 


t, i channel. 0-t only;l-t only;2-t and i 


seed 


Turn random number seed 



Table 2: Summary description of "param.dat" file. 



top quark FCNC interactions were introduced in MEtop through an effective la- 
grangian. Depending on which operators are "turned on", different physics will be gener- 
ated. Two different topologies are available: 2 — 1 — ?- 3 and 2 2 ^4.| The first one 
concerns "Direct top" production, and the second is related to "top+gluon" and "top+ 
light quark". 



Process Number 


Description 


Comments 


1 


Direct top (LO) 


Strong Op. only 


2 


top+gluon (LO) 


Strong Op. only; set PTmin 


3 


top+quark (LO) 


All Op.; set PTmin 


21 


Direct top (NLO) 


Strong Op. only 


22 


Direct top (NLO) + top+quark (LO) 


All Op. 



Table 3: Processes available in MEtop 



Strong FCNC top interactions are included in MEtop through two equivalent effective 
operators, one for the top and u-quark interaction, and the other for the interaction of 
the top with a c-quark. In process 1, only the strong coupling constants are needed. 

3 When "SpCorr" is set to 0, the top quark decay will not be performed in MEtop, that is, the 
generated events will have the topology 2 — > 1 and 2 — > 2. In this case the spin correlations are lost. 
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Process 2 has the same effective operators but due to the infrared divergences appearing 
in top+gluon production a cut in the top quark transverse momentum has to be set via 
the variable PTmin. In process 3, top + light quark @ production, all operators can 
contribute, strong, electroweak and 4F. It is now possible to choose which operators to 
include. Again a value for PTmin has to be chosen. 

Process 21 is inclusive direct top production at NLO and again only strong operators 
intervene. The NLO result is obtained by a matching procedure (as described previ- 
ously) which depends on one variable, PTmatch, to be chosen by the user. The cross 
section results are written in three files: "csDtopLO.txt", the LO result for direct top, 
"csDtopNLO.txt", the NLO increment relative to the LO result [o T N % - o T L °^ al ), and 
"csGtop.txt", the LO cross section for "top+gluon" process with a top quark transverse 
momentum above "PTmatch". Therefore the variable "PTmin" is irrelevant for this pro- 
cess. After the generation, the results are stored in one file in the Events folder named 
"DtopNLO.lhe" , containing 2 — > 1 — > 3 and the 2 — > 2 — > 4 configurations. These 
events constitute the inclusive direct top NLO event generation, and must subsequently 
be showered by PYTHIA using the PT-ordered scheme, in order to complete the match- 
ing procedure. Finally, with process 22, MEtop sums process 21 with process 3. The 
"PTmacth" variable plays the same role as in process 21 and "PTmin" will be the top 
transverse momentum cut, for the "top + quark" sub-processes. 

C.3 Running MEtop 

To run the package you just have to execute the command "./run.py" in the main direc- 
tory. Care should be taken when changing the values of the physical parameters and/or 
the process you wish to calculate. In such cases you must always recalculate the value 
of the cross section. In addition, if you change the process used for the generation, you 
must be sure that all cross sections pertaining the new process are calculated beforehand. 
This is mandatory because the generation is done using the cs*.txt files saved in the CS 
folder. 

C.4 Available Model files 

At the moment there are three different packages available in MEtop, with a different 
set of effective operators. The main reason to have the different packages is to make the 
generation of events faster. The lightest version includes only the strong sector. Then 
there are two other versions one with strong plus electroweak operators and the other one 
with strong plus 4F operators. 

In equation [32] we present the strong FCNC lagrangian as it is written in the package 
"MEtop_S_vxx.tar.gz" 

C s = col O uG + col O cG + h.c. (32) 

with 

O uG = i^u\ a a^(fu + hu l5 )tG% , O cG = i^c\W{fc + hc l5 )tG% (33) 
4 Here light quark stands for the set u, u, d, d, c, c, s, s, b, b. 
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and col, co2, fu, hu, fc, he are real constants to be chosen in the file param.dat. The 
constants f\, hi allow the choice of different chiralities while coj are overall normalization 
constants. Although it usually considered that ff + h? = 1, this relation has to be 
implemented by the user by a judicious choice of parameters and hi. 

The package "MEtop_SEW_vxx.tar.gz" contains the strong and electroweak sectors. 
The lagrangian introduced in this package is 

r — r j- ' ( rJi (nij i rJj cy-i _i_ ^Jo ryi i J^w) 

l-SEW - l-S + [ a uW U uW + a uB<f> U uB<t> + a <jm U <t>u + a <f>q U <t>q 

i,i=l,3 

+ (34) 
where the electroweak operators are 

Off = i^D^^im^r'qLj) , Ol^i(^D^)(u Ri ru Rj ) 

°uW = (<IU CV T I U Rj ) W lu , = («W <V M ^ ) 

and all coupling constants are real. In param.dat all coupling constant have the form 
coi. The relation between the coupling constants presented in equation [3H and the coj 
parameters to be chosen in param.dat is presented in table HI 



co3 ->■ atf w 


co4 -)• aft w 


co5 -> 


co6 ->■ a*^ 




co8 -> ag, 


co9 a£" tj 


colO 


11 (l,ut) 

coll-to^ 


co!2 c£« 


co!3 -> < 


col4 -> 



Table 4: Coefficient dictionary for Csew 
Finally, the file "MEtop_S4F_vxx.tar.gz" contains the strong and 4F sector 

£-SAF = £s + C\f u + (35) 

where the 4F lagrangians were presented in equations [29] and [301 The relation between the 
parameters in equations [29] and [30] and the corresponding coi parameters in the param.dat 
file is shown in table [5] 

Finally we note that any combination of parameters can be made in a new package 
and can be made available upon request. Generator and the different packages can be 
downloaded at http://coimbra.lip.pt/~miguelwon/MEtop/. 
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co27 a" 13 + a 11 , 13 

1 qq' 


co33 -> a" 23 + a 21 , 13 

' ^ qq 1 qq' 
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